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INTRODUCTION 


Aerospace structures are commonly fabricated using mechanically 
fastened joints. To better understand the complex behavior of such joints, it 
is important to develop accurate stress analysis methods. This is especially 
important for composite structures since they can be seriously weakened by 
fastener holes and often have rather complex failure modes. Within a multi- 
fastener structural joint, fastener holes may be subjected to the combined 
effects of bearing loads and loads that bypass the hole, as illustrated in 
figure 1. The ratio of bearing load to bypass load depends on the joint 
stiffness and configuration. 

The combined effects of bearing and bypass loads can be simulated by 
testing and analyzing single - fastener specimens (figure 1). However, very 
little research has been published on the stress analysis of single -fastener 
specimens subjected to combined bearing and bypass loading. In 1981, Ramkumar 
[1] and Soni [2] used a two-dimensional, finite element stress analysis to 
determine the stress state around the fastener in a single -fastener laminate 
subjected to bearing-bypass loads. The fastener load was modeled by imposing 
zero radial displacements on the load-carrying half of the fastener hole. In 
the same year, Garbo [3] used his BJSFM analysis to obtain the stresses in a 
single-fastener laminate subjected to bearing-bypass loads. The fastener load 
was simulated by specifying a radial stress boundary condition varying as a 
cosine over half the hole. 

In almost all practical applications, some clearance exists between 
the hole and the fastener. The presence of a clearance leads to a contact 
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region at the bolt-hole interface that varies nonlinearly with the load and 
influences the stress state around the hole. Although this nonlinearity 
complicates the analysis, it is important to make an accurate stress analysis 
to better understand the complex failure modes. Furthermore, for some 
combinations of compressive bearing-bypass loads the hole tends to close in on 
the fastener leading to "dual contact” . The effects of clearance have been 
presented by these authors in references 4-6. 

The objective of the present paper is to present a simple, direct 
s^ ess analysis method for a laminate with a clearance -fit fastener subjected 
to combined bearing and bypass loads in tension or compression, including 
dual contact, and then to study the effects of bearing-bypass loads on bolt- 
hole contact and local stresses. The present approach uses a linear- elas tic 
finite element analysis with an inverse formulation like that in Refs. 4-8. 
Conditions along the bolt-hole interface are specified by constraint 
equations that limit nodal displacements to a circular arc corresponding to 
the bolt diameter. These equations describe the contact conditions more 
realistically than the distributions usually assumed for radial displacement 
[1,2, 7, 8] or stress [3]. Furthermore, the present technique does not need an 
iterative- incremental method of solution which usually involves tedious node 
tracking along the contact arc [9-12]. The finite element analysis was 
performed using the MSC/NASTRAN computer code [13]. The material properties 
used in the analysis represent a quasi- isotropic T300/5208 graphite/epoxy 
laminate. The bolt was assumed to be rigid and the interface to be 
frictionless . 

Results are presented as curves relating the nominal bearing stress 
and the bolt-hole contact angle for various combinations of bearing and bypass 
loads. Also, the effect of bypass load on the bolt-hole contact angle is 
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presented for constant bearing loads. Hole boundary stresses for a range of 
bearing-bypass load proportions in both tension and compression are also 
presented for a typical bearing load level. Finally, a solution array for a 
range of bearing-bypass loads, typical load-contact variations, and local 
stresses are presented for dual-contact situations. 
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NOMENCLATURE 

radial clearance, m 

diametrical clearance, % 

hole diameter, m 

modulus in fiber direction, GPa 

modulus transverse to fiber, GPa 

shear modulus , GPa 

applied load, N 

bearing load, N 

bypass load, N 

polar coordinates, m, degrees 

bolt radius, m 

hole radius, m 

nominal bearing stress, MPa 

nominal net-section bypass stress, MPa 

specimen thickness, m 

x- and y-displacements , m 

specimen width, m 

Cartesian coordinates, m 
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bearing-bypass ratio, S^/S^ 
Poisson's ratio 

bolt-hole contact angle, degrees 
stress component , MPa 
radial stress component, MPa 
tangential stress component, MPa 


ANALYSIS 


The configuration and loading analyzed in the present study is shown 
in figure 2 for a typical tension bearing-bypass case. The gross applied load 
P is reacted partly in bearing P^ and partly as a bypass load P^ at the 
other end of the model. The nominal bearing stress and net-section 

bypass stress S are defined as shown in figure 2. The bolt was assumed to 
be rigid and the bolt-hole interface to be frictionless. The diametrical 
clearance c^ between the hole and the bolt was expressed as a percentage of 
the hole diameter d. The material represented in the analysis was a quasi - 
isotropic T300/5208 graphite/epoxy laminate with the following properties: 

E - 59.46 GPa and G - 20.4 GPa. Isoparametric, quadrilateral and triangular 
elements were used to model the laminate [4,13]. Nodes were placed at every 
0.9375 degrees along the hole boundary to model the contact angles accurately. 
Loads P and P^ were applied to the ends of the model. The desired 
bearing-bypass ratio f) was obtained by an appropriate choice of P and 

p b P * 

For a snug-fit joint (c^ - 0), the contact angle 6 ^ along the bolt- 
hole interface does not vary with and a simple linear stress analysis can 


4 



be used. However, for a clearance- fit joint (c^ > 0), the contact angle 8 ^ 
increases nonlinearly with as shown in figure 3. This nonlinear load- 

contact variation at the bolt-hole interface greatly complicates the stress 
analysis. Bypass loads also influence the nonlinear load-contact variation 
further complicating the analysis for combined bearing-bypass loading. 

Inverse Formulation 

The nonlinear load- contact variation shown in figure 3 could be 
accounted for in an iterative incremental scheme in which the nominal bearing 
stress is incremented in small steps and the corresponding contact angle 

8 ^ is determined iteratively at each load step. Such a procedure would 
require special purpose finite element programming and would also involve 
tedious node tracking along the contact arc [9-12]. Alternatively, an inverse 
technique [4] can be used in which a contact angle 8 ^ is assumed and the 
corresponding bearing stress is computed by a simple procedure which is 

described later. The inverse technique is simple to use because the boundary 
conditions for an assumed contact angle 8 ^ are fixed and known at the 
outset. For an assumed contact angle, the contact problem is linear, and the 
analysis procedure can be repeated for a range of contact angles to determine 
a nonlinear load-contact curve like that in figure 3. Therefore, although the 
contact problem is nonlinear, the inverse technique requires only linear 
finite element analyses; linear NASTRAN [13] procedures were used to solve 
this nonlinear problem. As previously mentioned, the conditions along the 
bolt-hole interface were specified by displacement contraint equations. The 
formulation of these constraint equations and the solution procedures are 
described in the following two sections. 
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Displacement Constraint Equations 


Figure 4 shows the analysis model for a typical case of contact at the 

bolt-hole interface. The bolt radius r, is smaller than the hole radius P 

b 

by the amount of radial clearance c. To simplify the analysis, the bolt is 
assumed to be fixed in space and the origin of the reference coordinate system 
is located at the center of the undeformed hole. Before load is applied, 
contact occurs only at point A. After the model is loaded on the two ends by 
P and , points along the hole boundary which lie within an assumed 

contact arc AB contact the frictionless surface of the fixed rigid bolt. 

Consider a point P(x,y) on the hole boundary within the assumed 
contact arc AB (see figure 4). Let u and v be the x- and y- 
displacements necessary to move point P from its original position to a 
point on the surface of the bolt. The deformed position of P may be 
described by the following equation: 


[(x - c) + u]‘ 


[y + v] = 


(l) 


By expressing x and y in polar coordinates, neglecting the higher order u 
and v terms, and noting that r^ - R - c, Eq. (1) may be rewritten as 
follows : 


A u + B v ~ C (2) 

where 

A * R cos 6 - c (3) 


B « R sin 6 (4) 

and 


6 



C - R c (cos 6 - 1) 


(5) 


Equation (2) is a constraint equation for the u- and v-displacements of any 
point P(x,y) on the contact arc AB. The quantities A, B, and C can be 
computed at the outset since they are functions of the initial geometry. In 
the finite element analysis, the displacements of each node within the contact 
region can be specified by applying Eq. (2) as a multi-point constraint [13], 
Note that the hole boundary beyond the contact arc is stress free. This fact 
will be used later in the analysis. 

As mentioned earlier, for some compression bearing-bypass loads, the 
hole tends to close on the bolt leading to dual contact at the bolt-hole 
interface. For example, point D on the hole boundary in figure 4 could move 
in the positive x-direction and make contact with the bolt surface. This 
would correspond to the onset of dual contact. Further loading would lead to 
dual contact as shown in figure 3. The contact angles 0 ^ and 9 ^ for dual 
contact would also, in general, vary nonlinear ly with load. However, the 
inverse technique can also be used for dual contact. Multi-point constraints, 
as described by Eq. (2), can be used to specify nodal displacements within the 
two contact regions. The corresponding combination of bearing and bypass 
loads can then be determined by the solution procedure described later. 

Solution Procedure for Single Contact 

The correct bearing stress S^ for an assumed contact angle 0^, was 
established using a simple procedure and the NASTRAN computer code. Point A 
(figure 4), which was assumed to be in contact with the rigid bolt, was fixed. 
Displacements of all nodes within the assumed^ contact arc were restricted to 
those allowed by the multi-point constraints given by Eq. (2). Thus all the 
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boundary conditions, including those along the bolt-hole interface, are 
defined. However, since the contact angle 9 ^ and the nominal bearing stress 
are nonl inear ly related, the correct corresponding to the assumed 9 ^ 

is still unknown. 

For a specified 9 the problem is linear and thus, the stresses in 

the plate are linearly related to and S . A linear equation relating 

o , S, , and S can be written as 
rr b np 

o rr (r,0) - F 1 (r ,9) S b + F^r.6) S np + F 3 (r ,9) (6) 


where F^, F 2 , and F^ are constants for a given r and 9. The first two 
terms in this equation represent the a components due to the applied 
bearing and bypass loads, respectively. The third term represents the a 
due to imposing the multi-point constraints and is a function of clearance, as 
indicated by Eqs . (2-5). 

For a given bearing-bypass ratio fi , we can write S as 


S - S, / /? 
np h ' y 


(7) 


Substituting Eq. (7) into Eq. (6) gives 


where 


a rr (r,0) - F 4 (r,tf) S b + F 3 (r ,9) 


F 4 (r ,9) « F 1 (r,tf) + F 2 (r 


( 8 ) 


(9) 


The hole boundary region beyond the contact angle 9 ^ is stress free. Thus, 
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d l < 6 < n 


( 10 ) 


* rr (R.*) 


This stress boundary condition was imposed at the end of the contact arc as 


a rr (R ’V “ 0 dU 

If is the correct bearing stress for an assumed contact angle 0^ t then 

from Eqs . (8) and (11), we have the following relation: 


F 4 (R ’V S b + - 0 


( 12 ) 


This equation can be rewritten as follows: 
S b “ -F 3 (R i « 1 ) / F 4 (R,^) 


(13) 


For a given 0, the values of F^R,^) and F^R,#^) were determined 

by the following procedure. For an assumed 6 the stress was selected 

arbitrarily, and a was calculated at the end of the contact arc using a 

finite element analysis. These and a values were then used in Eq. 

(8) to get one equation for F^R,^) and F 4 (R,0^). A second was 

selected (arbitrarily) and again the corresponding a at the end of the 

contact was calculated. The second set of S, and a values was used in 

b rr 

Eq. (8) to get a second equation for Fj(R t 0^) and F^(R,0^). The two 

equations were solved simultaneously to determine F^(R,0^) and F 4 (R,0^), 

which were then used in Eq. (13) to find S^. This procedure was repeated for 

ic 

a series of assumed 0 ^ values to determine the corresponding values. 

These pairs of 6 ^ and values can then be plotted to establish 
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versus 6 ^ curves for single contact in the clearance- fit joint. The above 
procedure was introduced and evaluated in reference 4. It can be successfully 
applied for both tension and compression bearing-bypass loads which lead to 
single contact. A slightly different procedure is used for dual contact and 
is described in the next section. 

Solution Procedure for Dual Contact 

As described earlier, compression bearing-bypass loads may cause dual 

contact at the bolt. The dual contact angles and 6 ^ (see figure 5) 

increase nonlinearly with load. However, for specified 6 ^ and the 

stresses in the plate are linearly related to S, , S . and the bolt-hole 

y b np 

contact stresses associated with the imposed multi-point constraints. A 
linear equation relating o S^, and S can be written for the end of 
each contact arc. For point B (see figure 5), 

a rr (R ^l } " K i< R ’V s b + V R ’ V S np + K 3 < R ’ V (14) 

and for point E, 

a rr (R ’*'V " s b + K 5 (R 9 w-» 2 ) S ^ + K 6 (R,tt -^) (15) 

The constants through are similar to those in Eq. (6). Again, the 

first two terms in both Eqs. (14) and (15) represent the a components 
due to the bearing and bypass loads, respectively. The third term represents 
the o associated with the imposed multi-point constraints. Equations (14) 
and (15) cannot be simplified further using Eq. (7), as was done earlier for 
the single-contact case (see Eqs. (6) and (8)), because the bearing-bypass 
ratio /9 that would lead to contact angles 6 ^ and 6 ^ is not known 
a priori. The constants through and through can be 
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determined by calculating and a rr ( R » n ~ ® ^° r t * iree different 

(arbitrary) combinations of and s n p* The calculated stresses and the 

corresponding values of and S^ p are then substituted into Eqs . (14) and 

(15) to yield three equations in K^, K 2 , and and three equations in K^, 

, and Kg. These six equations are solved simultaneously to determine 

through Kg. Now the correct values of and S corresponding to the 

assumed 6 ^ and # 2 values can be determined by imposing the boundary 

conditions at the ends of the contact arcs as a (R,0..) * 0 and o (R,7r-0„) 

rr v ’ 1 rr v 2' 

- 0. Substituting these boundary conditions into Eqs. (14) and (15) we have 
K 1 (R ’V S b + K 2 (R ^1 } S np + K 3 < R ’*i> - 0 (16) 

and 

K 4 (R,tt-S 2 ) S b + K 5 (R,w-* 2 ) S np + K 6 (R,7r-^ 2 ) - 0 (17) 

The only unknowns in the above equations are S, and S which can be 

d np 

determined by the simultaneous solution of Eqs. (16) and (17). The above 
procedure can be repeated for different sets of assumed contact angles 6 ^ 
and 9 ^ to determine the corresponding and S values. These results 

can be presented as a versus S plot as discussed in the following 

section. 

RESULTS AND DISCUSSION 

All of the results in this paper were obtained for the finite size 
plate shown in figure 2 with a clearance c^ of 1.2 percent. First, results 
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are presented for single contact with both tension and compression bearing- 
bypass loads. Next, the effect of bypass loading on bolt-hole contact is 
shown for a constant bearing loading. The corresponding hole boundary 
stresses are also presented. Results are then presented for dual contact. 


Contact angles are presented on a versus 

are presented for a dual -contact case. 


S 

np 


plot and local stresses 


Single Contact 

The effect of bearing-bypass ratio /J on the bearing stress-contact 
angle curves is shown in figure 6 for single contact. Note that f) - « 
corresponds to a tension bearing case with no bypass load and /? « -oo 
represents the corresponding compression bearing case. The bearing-bypass 
ratio /? was found to have a considerable effect on the bearing stress - 
contact angle behavior. For tension bearing-bypass at a bearing stress level 
of 400 MPa [6], a typical bearing strength for graphite/epoxy, the contact 
angle (figure 6) for /3 - 1 is about 30 percent larger than that for £ « « . 
Ir compression bearing-bypass at - 400 MPa, the contact angle for 0 - -1 
is about 25 percent smaller than that for - -«> . 

The variation of contact angle with bypass stress is shown in figure 7 

for three levels. Increasing the tensile bypass loading increased 6 ^ , 

while increasing the compressive bypass loading had the opposite effect. The 

small iog in the curve at S * 0 is caused by the small difference between 

np 

tension-reacted bearing and compression reacted bearing. For the - 400 
MPa case, dual contact initiated for a compressive bypass stress of about 
450 MPa. The secondary contact angle 9 ^ increased rather abruptly as the 
compressive S exceeded this value. Additionally, the decreasing trend for 
9 ^ reversed when dual contact developed. 
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For - 400 MPa, the tangential (a^) and radial ( £7 rr ) stresses 
around the hole boundary are shown in figure 8 for three different values of 
S for tension bearing-bypass loading. The peak value of the a ^ stress 


is not very sensitive to S 


np 


The peak usually occurs at 0 degrees for 


small values of S 


np 


However, the distribution of changes with 


increased S values, showing the increased contact angle as in figure 7. 
For S - 400 MPa, the a peak occurs at around 55 degrees. The increase 
in the peak a stresses rises nearly proportionally with the increase in 
S np stress. The location of the peak a ^ is usually a few degrees beyond 
the contact region. The location of the a .. peak is important because 
damage in composite joints often starts at this location [5]. 

The Ogg and hole boundary stresses for compressive 


np 


stresses are shown in figure 9. For S - 0 (f) - -®) , is mostly 


tensile. For S - -133 MPa, a 
np 
o 

(around 90 ) of the joint. For S 


np 


becomes compressive in the net-section 
- -400 MPa, the compression peak of a . 


is larger and more of the hole boundary is in compression. Unlike the tension 
bearing-bypass cases, the peak value of the a stress in compression 
bearing-bypass increases nearly proportionally with the bypass load and always 


occurs at 9 - 0 . For example the peak stress for S 


400 MPa is 


40 percent higher than the S 


np 


0 case. 


Dual Contact 

The S, and S values calculated for assumed values of 9 , and 
b np 1 

6 2 are plotted in figure 10. Each solid circular symbol represents a 
solution corresponding to assumed contact angles 9 ^ and 9^. The dashed 
curve represents the onset of dual contact. The region to the right of the 
dashed curve represents the dual-contact region. The double-dashed line 
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represents the p = -3 case. Note that this line runs almost parallel to the 
linear portion of the dual-contact threshold line (dashed line), indicating 
that the P - -3 case would never involve dual contact. Thus, dual contact 
would occur only when -3 < P < 0. The dashed-dot line represents the 
P - -0.5 case. It intersects the dual-contact curve (dashed) at S - -375 
MPa. Thus dual contact for this case would begin at S - -375 MPa 
(S^ = 187.5 MPa). Note that all the results in figure 10 were obtained for a 
clearance of 1.2 percent of hole diameter. A smaller clearance would cause 
dual contact to develop at a lower load level and vice versa. 

The effect of dual contact on the versus 9 ^ curve is shown in 

figure 11. The versus 9^ curve before dual contact was generated using 

the procedure for single -contact analysis described earlier. The remainder of 
the curve was constructed from figure 10 by selecting S^, 9^ t and 9^ values 

on the P - -0.5 line. Some interpolation was used. The versus 9 ^ 

relationship changes after the onset of dual contact. The versus 9 ^ 

curve is also nonlinear and 9 ^ increases more rapidly than 9 ^ with 
increasing load. 

An important consequence of dual contact is that it allows load 
transfer across the bolt and, therefore, reduces the stress concentration 
around the fastener hole. Figure 12 illustrates this effect by comparing 
hole boundary stresses for a case with dual contact to stresses for an open 
hole case at the same load level. The peak o ^ for the open hole case is 13 
percent higher than that for the dual-contact case. This suggests that dual 
contact is advantageous since it reduces stresses around the hole and 
increases the joint strength. Smaller bolt-hole clearances promote dual 
contact and, therefore, should produce higher joint strengths in compression. 
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CONCLUDING REMARKS 


A simple method has been developed for the stress analysis of a 

laminate with a clearance -fit fastener subjected to combined bearing and 

bypass loading in tension or compression, including dual-contact situations. 

This method uses a linear elastic finite element analysis with an inverse 

formulation. The present method is simple to apply and can be implemented 

with most general purpose finite element programs. The method was applied to 

study the effects of bearing-bypass load proportioning on the bolt-hole 

contact angles and local stresses. The initial clearance between the smooth, 

rigid bolt and the hole was 1.2 percent of the hole diameter in all analyses. 

Material properties for the plate represented a quasi -isotropic, 

graphite/epoxy laminate. The bearing-bypass proportions were expressed in 

terms of fi , the ratio of bearing stress to bypass stress S . 

The bearing-bypass ratio fi was found to have a considerable effect 

on the nonlinear contact behavior. Under tension bearing-bypass with 

- 400 MPa, the contact angle was 30 percent larger for fi - 1 than that for 

fi - ». Compressive bypass loading had the opposite effect, the contact angle 

was 25 percent smaller for fi - -1 than that for fi - -<*> at S, - 400 MPa. 

b 

For single -contact situations under tension bearing -bypass , the peak 

tangential stresses around the hole boundary increased proportionately with 

the S stress. The peak radial stress was not very sensitive to S 

np r J np 

However, under compressive bearing-bypass loading, both the peak tangential 
and peak radial stresses were considerably influenced by S . For P 
greater than -3, the hole tended to close on the fastener leading to dual 
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contact. Dual contact allows load transfer across the fastener and, 
therefore, reduces the stress concentration around the hole. Dual contact, 
therefore, has a beneficial effect on the joint strength. 

These results illustrate the general importance of accounting for 
bolt-hole clearance and contact to accurately compute local bolt-hole stresses 
for combined bearing and bypass loading. 
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Figure 1. - Bearing-bypass loading within a multi- fastener joint. 
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Figure 2. - Plate configuration for bearing-bypass loading. 
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Figure 3. - Nonlinear relationship between bearing stress and contact 

angle . 
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Hole clearance notation for single-sided contact 
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Figure 5. - Contact angle notation for dual 
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Figure 7. - Variation of contact angle with bypass stress. 
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Figure 8. - Stresses along hole boundary for tension loading. 
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Figure 9. - Stresses along hole boundary for compression loading 
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Figure 10. - Solution array for dual contact. 
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pression. The results illustrate the general importance of accounting for bolt-hole clearance and and contact 
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